Calibrating dipolar interaction in an atomic condensate 
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We revisit the topic of a dipolar condensate with the recently derived more rigorous pseudo- 
potential for dipole-dipole interaction [A. Derevianko, Phys. Rev. A 67, 033607 (2003)]. Based on 
the highly successful variational technique, we find that all dipolar effects estimated before (using 
the bare dipole-dipole interaction) become significantly larger, i.e. are amplified by the new velocity- 
dependent pseudo-potential, especially in the limit of large or small trap aspect ratios. This result 
points to a promising prospect for detecting dipolar effects inside an atomic condensate. 
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Interactions make life interesting. To a large degree, 
they determine both the kinematic and dynamic prop- 
erties of a physical system. In recent years, atomic 
quantum gases have become testing grounds for inves- 
tigating interaction effects. At typical temperatures for 
these quantum gases, the dominant interaction, the bi- 
nary atomic elastic collision, is isotropic and thus char- 
acterized by the s-wave scattering length a ss . The ma- 
nipulation of its strength from strong to weak, and its 
character from attractive to negative with a Feshbach 
resonance, has become one of the continuing highlights. 
At a more detailed level, however, atoms are composite 
particles, e.g. possessing magnetic dipole moments, from 
the electron and the nuclear spin. The resulting dipo- 
lar interaction between atoms, is anisotropic, and consti- 
tutes an exciting new development. While much weaker 
as compared to typical isotropic s-wave interaction, its 
experimental detection is only a matter of time, consid- 
ering the rapid pace of advances in this active field. 

For a condensate of N atoms interacting via a potential 



polarized atoms (along z) including dipolar interaction 



the total energy functional is 
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where ip{r) is the condensate wave function and R = 
f- f'. V cxt (r) = Mw 2 p {x 2 + y 2 + X 2 z 2 )/2 is the trap 
potential, assumed harmonic and of axial symmetry with 
radial (axial) trap frequency uj p (lu z — Xlu p ). 

The real (bare) potential V(R) in the interaction en- 
ergy [the 2nd line of Eq. |Q] is usually replaced by a 
pseudo-potential V, which for an isotropic short ranged 
interaction takes the contact form 



V(R)=g6(R), 



(2) 



with g = 4irh 2 a ss /M . a ss is the s wave scattering length 
of V(R). To date, this pseudo-potential approach has 
proven remarkably effective for most studies. 

In this Letter we revisit the topic of a condensate of 



Vub{R) = ^§P 2 (cos0 fl ), 
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where Or is the polar angle of R, and P2O) is the 2nd or- 
der Legendre polynomial. This problem is important be- 
cause the non-spherically symmetric interaction Eq. © 
leads to interesting low energy collisions due to the pres- 
ence of both the 'short-' and 'long-' range characters 
We have previously suggested a pseudo-potential 
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by matching its first Born amplitudes to the complete 
scattering amplitudes of V(R) [2| and confirmed its va- 
lidity, provided the dipole moment is not much larger 
than a Bohr magneton (/J-r) and the collision is away 
from any shape resonances 

Atomic dipolar condensate with interaction J3J has 
been studied by many groups. A lot has been learnt 
about its ground state and the associated stability 0, 
HHI1], collective excitations 0, 0, H, , free expansion 
dynamics 0, , and the potential existence of several 
exotic phases in an optical lattice [Lq|. Recently, it was 
discovered that the ground state density profile remains 
an inverted parabola in the Thomas- Fermi limit fl3| . 

Following Huang and Yang for spherically sym- 
metric potentials, Derevianko recently proposed a more 
rigorous pseudo-potential V ans i applicable to anisotropic 
interactions including regions near collision resonances 
Due to its dependence on the relative momentum, 
the interaction energy of V ans i is most conveniently ex- 
pressed in momentum representation as [TH H^ ] 



4 / dkxdkdk'4>*(ki)(t>*(ki - 2k)v(k, k') 
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where 4>(k) is the Fourier transform of ip(r). k and k ' are, 
respectively, one half the pre- and post-collision relative 
momenta of the colliding pair. The momentum represen- 
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tation of V ans i, denoted by v(k, k') takes the form If 
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with a s( j the generalized scattering length due to the cou- 
pling of the s (I = 0) and d (I = 2) partial wave chan- 
nels. Both a ss and a s( j are obtained from the zero energy 
T-matrix elements of a coupled multi-channel scattering 
calculation pj. In the low energy limit [l5j . 

T DD {k,k') = VEP 2 (cos6 k >) + P 2 (cos0 fe ), 

where the 1st term is momentum-independent, essen- 
tially corresponds to the bare dipolar potential Eq. 
the 2nd term, on the other hand, depends on the mo- 
menta. It arises from the rigorous construction of the 
pseudo-potential. The main purpose of this study is to 
calibrate how it modifies the properties of a dipolar con- 
densate. To our surprise, we find previous studies with- 
out the momentum-dependent 2nd term have severely 
under-estimated the strength of the dipolar interaction. 

Although Vansi is non-Hermitian, a variational study 
can nevertheless be performed within the mean field the- 
ory Vansi is constructed by matching its (two-body) scat- 
tering solution in the asymptotical (R large) limit to that 
of the real potential V(R) ^5|. It is non-Hermitian be- 
cause its scattering solution differs from the real one in 
the short-range. Within the mean field approximation, 
however, all condensed atoms share the same spatial or- 
bital, a smoothed or "coarse grained" condensate wave 
function (valid over length scales much larger than the 
range of atomic interactions) . Thus the mean field theory 
is limited to a subspace of nonsingular functions of the 
complete Hilbert space, where the pseudo-potential V an si 
is Hermitian and leads to an unitary time evolution. As 
expected, with a Gaussian ansatz 
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of variation parameters w„ and w z , we find 
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indeed being real (see also Ref. k = w p /w z is the 

condensate aspect ratio and x{ K ) — Xo{ K ) + Xi{ K ) with 



Xo(k) 
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1 - 3k 2 H{k) 



-( k 2 -1)H(k), 



and H{k) = tanh" 1 y/l - k 2 /VT 



The Xo( K ) term 



arises from the bare dipolar interaction Eq. @ and is 
known before [Hill, while the xi( K ) term is due to the 
momentum-dependent second term of TDD{k,k'). 
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FIG. 1: The function £(«). 



Both Xo( K ) an d Xi( K ) are nionotonically increasing 
functions of k and vanish at k = 1. More specifically, 
Xo is bounded between —1 and 2, while \i diverges to 
— oo and oo at k — and oo. Thus the net dipolar ef- 
fects from the pseudo-potential Eq. © are larger, or 
more prominent than realized before using Vpp" (R) . As 
a comparison we plot in Fig. = Xi{ K )/Xo( K )- We 

note is relatively flat (~ 4) for a condensate with 
moderate aspect ratio k ~ 1, i.e. the pseudo-potential 
Eq. © differs from Eq. only by a scaling factor. 
Therefore the pseudo-potential V^ n (R) remains valid 
near k ~ 1, provided C3 is proportionally scaled (renor- 
malized) by a factor of ~ 5. This enhancement clearly 
points to the effect of off-shell k ^ k' collisions. For on- 
shell collisions [k = k' in Tod (k, k')] , approximately a 
factor of 2 enhancement arises due to the symmetriza- 
tion of the scattering amplitude. In the extreme limits 
of k, collisions are restricted to either ID or 2D where 
different scattering behavior may arise 0]> making the 
use of Eq. © questionable because a s dis obtained from 
the zero energy (k — ► 0) 3D T-matrix [jj . 

Using a p = yfh/Mbjp (fkjp) as unit for length (energy), 
the dimensionless form of energy per atom becomes 
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where q p _ z — w PiZ /a p . The contact and dipolar inter- 
action parameters are V ss = W 71 " Na ss /a p and V s d — 
— \/5/(27r) Na s d/a p . The condensate widths (q p o,q z o) 
are obtained through a minimization according to 
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FIG. 2: The stability diagram of a dipolar condensate when 
Pss = 0. 



10 




(c) 



5 \<* 



"2 







10 



10 



2 V 



FIG. 3: The stability diagram for a dipolar condensate with 
a ss > 0. (a) and (b): the A dependence of for V ss = 1 
(solid line), 10 (dashed line), and oo (dash-dotted line); (c) 
and (d): the A dependence of P s ( s +) (c) and Pi" 5 (d) for C = 
±0.5 (solid line), ±1 (dashed line), and ±2 (dash-dotted line). 



where kq = 1po/<]zO and 
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A solution (q p o, q z o) is stable if 
~8 2 £ d 2 £ 



dql dq 2 z 



d 2 £ 



> 0. 



(13) 



We note A is tunable as mentioned before, one can 
also vary a ss with a Feshbach resonance [2(J, |^| • Close 
to dipolar induced shape resonances, a s d can be similarly 
tuned to large or small and positive or negative 1] . The 
positive valued a s( j can be easily understood by consid- 
ering the interaction between two polarized dipoles. If 




FIG. 4: (a) The dependence of fe s t on a BB for A = 6.8/17.35. 
The dashed line is the result without dipolar interaction, (b) 
The dependence of fe s t on A for a ss — — 0.5(as) (dashed line), 
— Has) (dash-dotted line), and — 4(as) (dotted line). The 
solid line is the result without dipolar interaction. 



they were placed in a plane perpendicular to their po- 
larization (tt), they attract each other; while they repel 
each other when placed along the direction of their po- 
larization (—>—»•), (note the difference with respect to the 
bare dipolar interaction 0). Based on these consider- 
ations, we proceed to study the properties of a dipolar 
condensate by varying V ss , V s a, and A. 

Figure [21 shows the stability diagram of a dipolar con- 
densate when a ss = 0. For a s( j < 0, the condensate is 



stable if V sd < V. 



(+) 

"sd 



and there exists an 'always sta- 
7.56, which is greater than the 
Interestingly, for a sc j > 0, the 
>(-) 



ble' region if A > A 
previous estimate 

condensate is stable if V s & > V^' and it is always stable 
for A < A_ ~ 0.142. 

When a ss > 0, it is convenient to define C = 'P s d/'P SS 
which measures the relative strength of the dipolar in- 
teraction. C can be changed by tuning a s( j. As shown in 
Fig. 01(a) and (b), for a given a ss > 0, the condensate 
is stable if C*- - ' < C < C^ + K Similar to the previous 
result 0, the critical value of A for the 'always stable' 
region is approximately C-independent and equal to A± 
for a ss — 0. For a fixed C, the condensate stability can 
be modified by tuning V ss - From Fig. El(c) and (d), we 
see that only when V ss < for C > and V ss < 
for C < is the condensate stable. Varying V ss for a 
constant C can be achieved by changing AT. 

When a ss < 0, the condensate becomes unstable if 
N exceeds a critical value iV cr even without dipolar 
interaction. This instability is typically measured by 
the stability coefficient k st = N cr \a sc \ / a^o p2l |23| with 
dho = y/h/muj and Co = (a^ovu^) 1 / 3 . To study the dipo- 
lar induced modification of fc s t, we consider the experi- 
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FIG. 5: The N dependence of A fc for a BB = 5(a B ) (a), 0.5 (b), (c), -0.5 (d), -5 (e), and A = 0.2 (solid line), 0.6 (dashed 
line), 2 (dash-dotted line), 6 (dotted line). 



ment jH E3| of 85 Rb atoms in state \F = 2,M F = 2) 
with a magnetic dipole moment of /i = 2^b/3. We 
take a s d as obtained for the multi-channel scattering cal- 
culation PJ. For Up — (2tt) 17.35 (Hz), we then have 
V sd « 5.0 x 10~ 6 N. Figure H (a) shows the a ss de- 
pendence of fc s t for A = 6.8/17.35. As before dipolar 
interaction destabilizes a condensate for this configura- 
tion; However, the enhanced effect due to the more rigor- 
ous pseudo-potential Eq. ijBJ is more profound: even for 
a ss < —3 (as), the effect of dipolar interaction is still visi- 
ble. In Fig. 0](b), w e plot the A dependence of k st . Since 
"Psd > 0, dipolar interaction destabilizes the condensate 
at small values of A; when A is large, dipolar interaction 
can also stabilize an otherwise attractive condensate. 

Finally, we briefly consider the condensate aspect ratio 
kq. We define its relative change |6j 

Ak = *o(7W0)- "0(7^ = 0) (w) 

K (Vsd =0) 

as a measure of whether it is possible to detect the dipolar 
interaction from imaging condensate shape within cur- 
rent experiments. Figure [S] shows its iV-dependence for 
various a ss and A. For certain A, A K can be as high as 



20% even with N = 2 x 10 5 , if |a ss | is tuned small. 

In conclusion, we have calibrated the properties of a 
dipolar condensate using the more rigorous anisotropic 
pseudo-potential [15j, based on a variational calculation. 
Significant enhancement were found to all dipolar effects 
predicted previously. In the limit of weak s-wave inter- 
actions, when P ss < 1 due to a small a ss and/or a small 
N, our results are clearly valid based on previous stud- 
ies with the same variational technique 0, H, EH D3 03 , 
where extensively comparisons were preformed to justify 
the trial function Eq. J7J) fl7L ITsj . including the pres- 
ence of a weak dipolar interaction 0, 0, ^J. We note 
this is also the interesting limit where experimental de- 
tection of coherent dipolar interactions will likely occur. 
The variation approach gives reliable stationary conden- 
sate properties, although not the density profile itself, 
even when the interaction energy is large [lg. The re- 
sults from P ss > 1, although consistent, needs further 
improvements with more accurate numerical calculations 
(over the present variational method). 

We acknowledge insightful discussions and communi- 
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NSF and NSFC. 
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